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We study nonlinear inter-subband microwave absorption of electrons bound to the liquid helium 
surface. Already for a comparatively low radiation intensity, resonant absorption due to transitions 
between the two lowest subbands is accompanied by electron overheating. The overheating results 
in a significant population of higher subbands. The Coulomb interaction between electrons causes 
a shift of the resonant frequency, which depends on the population of the excited states and thus 
on the electron temperature T e . The latter is determined experimentally from the electron pho- 
toconductivity. The experimentally established relationship between the frequency shift and T e is 
in reasonable agreement with the theory. The dependence of the shift on the radiation intensity 
introduces nonlinearity into the rate of the inter-subband absorption resulting in bistability and hys- 
teresis of the resonant response. The hysteresis of the response explains the behavior in the regime 
of frequency modulation, which we observe for electrons on liquid 3 He and which was previously 
seen for electrons on liquid 4 He. 

PACS numbers: 67.90.+Z, 73.20.-r, 78.70.Gq, 42.65.Pc 
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I. INTRODUCTION 

Free electrons outside liquid helium can be trapped at 
the vapor-liquid interface owing to the attractive polar- 
ization potential above the surface and a repulsive barrier 
at the surface, which prevents electrons from entering 
the liquid^— As a result, the electron motion perpen- 
dicular to the helium surface is quantized, with energies 
e a = —R/a 2 (a— 1,2,..), where R is the effective Rydberg 
energy, R w 7.6 K for liquid 4 He and R ~ 4.2 K for liquid 
3 Hc. An external electric field E± applied perpendicular 
to the surface Stark-shifts the energy levels. The shift 
depends on the quantum number a. Consequently, the 
level spacing can be varied by changing E±. If there are 
no other confining fields the electrons are free to move 
along the surface and their overall energy spectrum con- 
sists of two-dimensional (2D) subbands. 

Resonant inter-subband optical absorption due to tran- 
sitions between the subbands has been directly observed 
for surface electrons (SEs) above liquid 4 He using ra- 
diation in the millimeter— £ and far-infrared^ ranges. 
Recently resonant absorption of millimeter-range mi- 
crowaves was also observed for electrons on liquid 3 HeZ 
In the experiments, SEs were tuned in resonance with 
radiation by adjusting the value of E± and optical 
absorption was detected as the variation in radiation 
power passing through the system. Inter-subband ab- 
sorption has been also probed indirectly by observing a 
change in the conductivity of irradiated SEs on both 4 He 
and 3 Hei&2 This change is mainly caused by radiation- 
induced overheating of the electron system.— The over- 
heating occurs both where electrons are scattered pri- 
marily by helium vapor atoms or by surface vibrational 



modes, ripplonsi^i — 

The recent interest in resonant absorption was stimu- 
lated by the proposal to use laterally confined electrons 
on the helium surface as qubits for quantum computing, 
with the qubit states being the out-of-plane states with 
a = 1 and 2 i 15 i 16 Resonant interaction with an exter- 
nally applied microwave field provides a natural way to 
manipulate the qubit state. Coupling between the qubits 
comes from the Coulomb interaction between the elec- 
trons. The interaction energy depends on the electron 
quantum numbers a, because the electron distance from 
the surface depends on a. Therefore, the change of the 
quantum state of one electron results in the shifts of the 
quantum levels of other electrons. 

In microwave absorption experiments with electrons 
that are free to move along the surface, one would ex- 
pect to see a shift in the resonance frequency for the 
transition from the ground to the first excited sub- 
band (and also in other transition frequencies) depending 
on the population of excited subbands. The subband- 
population dependent shifts of the energy levels due to 
many-electron effects were predicted earlier for semicon- 
ductor heterostructuresi 17 ' 18 In such systems the electron 
correlation effects are a perturbation, because the system 
is close to metallic. 

A distinctive feature of the electron system on he- 
lium is that it is strongly correlatedi 2 i 3 i 19 ~ — Therefore 
the frequency shifts turn out to be much stronger than 
in semiconductors. Signs of the dependence of the fre- 
quency shift on the state of the system were seen earlier 
in a complicated response to radiation with modulated 
frequency 2 ^ and the power-dependent shift of the photo- 
conductivity signal^ However, no systematic studies of 
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these effects have been carried out. 

The state-dependent shift of resonance frequency in- 
troduces nonlincarity in the absorption rate. Such non- 
linearity can cause bistability and hysteresis of resonant 
response. Hystcretic effects related to nonlinear absorp- 
tion are familiar already from the early work on ferro- 
magnetic resonance in ferrite crystal disks^~— For 2D 
systems, magnetization bistability was observed in spin- 
polarized atomic hydrogen gas absorbed on the supcrfluid 
helium filmj^l. For semiconductor heterostructures, the 
nonlincarity of inter-subband absorption was predicted 
and studied since late 80s2£~— and absorption bistability 
has been long sought 

In this work, we present theoretical and experimen- 
tal studies of nonlinear absorption effects for a strongly 
correlated electron system on liquid helium surface. A 
short account of the results appeared earlier^ We pro- 
vide a theory of the frequency shift associated with elec- 
tron correlations and show that this shift causes bista- 
bility and hysteresis of resonant response. We describe 
experimental results on the frequency shift and nonlinear 
absorption effects for electrons on liquid 3 He in the tem- 
perature range 0.2—0.4 K that covers both vapor-atom 
and ripplon scattering regimes. The response is observed 
as a change of the longitudinal conductivity <j xx of SEs 
in a weak perpendicular magnetic fields. The frequency 
shift is measured as a function of electron temperature 
T e , which is determined from the relative change of a xx 
at resonance. The shift is compared with the theoretical 
estimate and a reasonable agreement is found. 

At high values of input microwave power, the absorp- 
tion cross-section shows a jump as a function of E±_ and 
depends on the direction of sweeping Ex, that is, dis- 
plays hysteresis. This behavior agrees well with our the- 
oretically predicted bistability of nonlinear response. We 
also describe the results of the experiment, in which res- 
onance is observed by direct measurements of the mi- 
crowave power absorbed by the electron system. We show 
that the complicated behavior of the response to radia- 
tion with slowly modulated frequency can be explained 
by taking into account the nonlinearity of the absorption 
rate. 

In Sec. II we derive a many-body quantum kinetic 
equation for a strongly correlated system and use it to de- 
scribe nonlinear resonant absorption and to predict bista- 
bility of resonant response. In Sec. Ill we describe the 
experimental techniques. In Sec. IV the experimental re- 
sults are presented. In Section V we discuss the results 
of the measurements and make comparison with the the- 
ory. The experimentally observed nonlinear effects, in 
particular the bistability and hysteresis of inter-subband 
absorption, are also discussed in this section. Section VI 
provides a summary of the main results and the conclu- 
sions. 



II. THEORETICAL BACKGROUND 

A. Qualitative picture of electron dynamics 

In the neglect of the electron-electron interaction, elec- 
tron motion normal to the helium surface is described by 
the Schrodingcr equation 

h 2 d 2 

- -^-^2^ a (z) + V{z)ip a (z) = e a Tp a (z), (1) 

where z is the out-of-planc coordinate, ip a (z) and e a 
are the eigenfunctions and eigenvalues (a = 1,2,...), 
m is the electron mass, and V(z) is the electron po- 
tential energy, V(z) = —A/z + eE±z for z > 0. Here, 
A = e 2 (e — l)/4(e + 1), — e is the electron charge, e is the 
dielectric constant of liquid helium, and Ex is the elec- 
tric field applied normal to the surface. On the surface, 
z = 0, the potential V(z) has a high step ~ 1 eV. Neglect- 
ing penetration of the electron wave function into the liq- 
uid, Eq. (fTJ) can be solved numerically with the boundary 
condition ip a (0) = 0. For the field Ex < 10 2 V/cm, the 
typical localization length of the low-lying electron states 
is r B = h 2 /Am w 0.8 x 10" 6 cm" 1 for 4 Hc. 

Of central interest for this paper is electron response 
to a microwave field which is resonant with the a = 1 — >• 
a = 2 transition. This response depends on electron re- 
laxation. To describe the relaxation we note that, for 
electron densities n s and temperatures T studied in ex- 
periment, the electron-electron interaction is strong. The 
interaction is characterized by the plasma parameter Y 
and the characteristic plasma frequency cj p , 

Y = e 2 (7m s ) 1 / 2 /k B T, lu p = (2 7 re 2 n 3 / 2 /m) 1 / 2 . (2) 

Most of the existing experimental results on electrons on 
helium refer to the range r>l. This means that the 
electron system is strongly correlated, SEs form a non- 
degenerate liquid or, for Y > 130, a Wigner solid; the 
experiment described in this paper refers to the range 
where the system is a liquid. The reciprocal plasma fre- 
quency w" 1 gives the characteristic time of inter-electron 
in-plane momentum and energy exchange. In our exper- 
iment uip > 6 x 10 9 s^ 1 . This was the shortest relaxation 
time in the range of temperatures studied in our experi- 
ment. Also, uj p exceeded the microwave absorption rate, 
see below. 

The fast inter-electron energy exchange leads to a ther- 
mal distribution over the in-plane electron energy, with 
an effective electron temperature T e . We note that the 
time it takes to form such distribution in a correlated 
electron system can be longer than w" 1 . For example, 
if one thinks of electrons forming a Wigner crystal, this 
is the thermalization time of the phonons of the Wigner 
crystal. Still this time is much shorter than the char- 
acteristic energy relaxation time due to inelastic electron 
scattering by helium vapor atoms or ripplons/phonons in 
helium. Of major interest to us is the parameter range 
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where k^T e > ftio p . In this range electron in-plane mo- 
tion is semi-classical and the electron in-plane state can 
be described by well-defined momentum and coordinate. 

The temperature T e should be the same in all sub- 
bands, even though the distribution over the subbands 
does not necessarily has to be thermal. This is a con- 
sequence of the strong difference between the character- 
istic in-plane and out-of-plane inter-electron distance in 

]/2 . 1/2 

the correlated system, n s and rg- Since n s 3> rg, 
the in-plane momentum exchange Ap ee is essentially the 
same whether the electrons are in the same or in different 
subbands, leading to the same in-plane momentum distri- 
bution irrespective of the subband. At the same time, the 
electron-electron interaction practically does not lead to 
inter-subband transitions, because it requires short-range 
collisions which do not occur in the correlated system. 

Other scattering mechanisms come from the interac- 
tion of SEs with the surrounding, primarily with he- 
lium vapor atoms, for comparatively high temperatures, 
T > 0.3 K for electrons on 3 He, and with surface 
waves and bulk excitations in liquid helium, for lower 
temperatures^^ Scattering by vapor atoms and by sur- 
face capillary waves, ripplons, is mostly quasi-elastic. 
This is a consequence, respectively, of the large ratio of 
the helium atom mass to the electron mass and the slow- 
ness of ripplons: even for the ripplon wave number as 
large as 1 the ripplon frequency is very small, ~ 10 s Hz, 
and the ripplon energy is much smaller than the average 
energy of in-plane electron motion. 

B. Hamiltonian of the many-electron system 

— 1/2 

Using the condition n s 3> TBi the Hamiltonian of 
the isolated electron system can be written as Ho + Hel\ 
where the leading term is 

H = H K + H<® + H z , JJ K = iV p 2 n , 

2m L — ' 

n 

#ie 0) = yl>nn" H z = J2^ a - (3) 

Here, r„ and p n ate lateral coordinate and momentum 
of an nth electron, respectively, r nn i = |r n — r n /| is the 
in-plane inter-electron distance, and 

af = \a) nn {f3\, (4) 

where \a) n = [ip a ]n is the ath out-of-plane state of the 
nth electron. 

The part of the electron-electron interaction incorpo- 
rated into Hq corresponds to the approximation of the 
electrons being in a plane. The part of the interaction 
related to the out-of-plane electron displacement is de- 
scribed by HeV . To the leading order in r-Q/r nn i 

H ee=-\ E(^-*"')>L'- ( 5 ) 



As we will see, this interaction leads to an ana- 
log of the depolarization effect in semiconductor 
heterostructires . 1 7 i 18 1 29 ~ ^ However, because SEs on he- 
lium are strongly correlated the effect is much stronger 
and is described in a qualitatively different way. 

The electron coupling to a resonant microwave field is 
described by the Hamiltonian 

H F = -hn R cosoj F t K 1 / 2 + al 1 ) . (6) 

n 

Here, fin = eEuw^/h is the Rabi frequency; Emw 
and uip are the amplitude and frequency of the radiation 
field, respectively, and z a p = (a|z|/3). Frequency top is 
assumed to be close to the transition frequency of the 
electron system, CJ21 = (e2 — e i)/fij l w 2i — w_f| -C CJ21- 

The nonlinear response of a strongly correlated system 
can be analyzed using the many-electron density matrix. 
It depends on the states of in-plane and out-of-plane mo- 
tion of all electrons. It is convenient to write it in the in- 
teraction representation using canonical transformation 
U(t) = exp(-iHtt), where 

H t = H + h(uo F - w 21 ) J2 <ff + H He ■ (7) 

n 

The term Hne in Eq. ([7|) describes excitations in liquid 
helium and the almost ideal gas formed by helium vapor 
atoms. The transformed Hamiltonian of coupling to the 
field U'HpU is independent of time in the rotating wave 
approximation. We emphasize that the electron-electron 
interaction H^) is incorporated into H tl this is not a 
perturbation. 

We write the density matrix in the interaction rep- 
resentation with Hamiltonian H t using the Wigner 
representation with respect to in-plane motion as 
Po(t; {a„, P ni v n , p n })- For weak coupling to ripplons and 
helium vapor atoms, one can obtain a many-electron ki- 
netic equation for po- This equation is Markovian in 
time slow compared to /i/fcsT e , o;^ 1 and presents an im- 
mediate extension to the multi-subband case of the ki- 
netic equation discussed previously in the one-subband 
appr oximat ion . 3 9 

Because the Coulomb interaction is strong, to zeroth 
order in the interaction with excitations in helium and 
in H^) function po depends on the coordinates and mo- 
menta of individual electrons in terms of the total in- 
plane momentum P = ^ n p„ and the total in-plane en- 
ergy Hk + H^e ■ The distribution over the total in-plane 
energy is of the Boltzmann form. Moreover, in the classi- 
cal range k^Te > hu p the distribution over momenta p„ 
is also of the Maxwell form. 

Coupling to helium vapor atoms and ripplons causes 
mixing of states with different P and also of different 
out-of-plane states \a) n . The coupling Hamiltonian can 
be written as 

ff, = EE^ e ' ,r "<' ( g ) 

q,n a, /3 
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Here, V is an operator with respect to the variables of 
ripplons or phonons in helium and the positions of the 
helium vapor atoms, while Hq is the lateral momentum 
transferred from these excitations to an electron. Terms 
with a = ft describe intrasubband coupling, whereas 
terms with a ^ (5 describe mixing of states in different 
subbands. 

The characteristic values of the transferred momentum 

— 1/2 

are comparatively large, q 3> n s , because the density 
of states of excitations in helium and often the interac- 
tion strength increase with q. Thus scattering by helium 
excitations is short-range. As a result, in the electron 
liquid each electron is scattered individually, processes 
where two electrons are scattered by the same ripplon, 
phonon, or a vapor atom can be disregarded] 19 ' 39 Since 
between the scattering events the electrons have time to 
exchange lateral momentum and energy with each other, 
the scattering rates can be averaged over lateral electron 
motion with the Boltzmann factor. 

Resonant response of the electron system to microwave 
field is independent of the electron lateral motion. Of 
interest for the study of this response is the single-particle 
density matrix 



Pa/3 



(t) = Tr [a^p (t)] = AT 1 £ Tr [a^p (t)] . (9) 



The trace is taken over the in-plane coordinates and mo- 
menta and over the out-of-plane states of all electrons, N 
is the total number of electrons. In a spatially uniform 
system p a p is independent of the subscript n in Eq. ([9]). 



C. Many-electron shift of the transition frequency 

Qualitatively, the major effect of the interaction h2^ 
on resonant absorption is the linear Stark shift of the 
frequency of the inter-subband |1) — > 2) transition of an 
electron due to the out-of-plane component of the elec- 
tric field created by other electrons. This field for an 
nth electron is given by e~ 1 dH < ^) / dz n . As seen from 
Eq. (|S|), the field depends on the relative distance be- 
tween the electrons normal to the surface and therefore 
on the states of out-of-plane motion of different electrons, 
as well as on the inter-electron distance in the plane. We 
emphasize that we are interested in the field on an elec- 
tron, not in the average field in the electron layer, which 
is significantly different. 

Formally, the term H^) couples a single-particle den- 
sity matrix, in particular that given by Eq. ([9]) , to a two- 
particle one. This coupling is described by the term in the 

kinetic equation (dtf>o)ee = ift -1 Poit), flee (i) , where 

Hil\t) = W(t)H$U(t). The density matrix in the in- 
teraction representation slowly varies in time. Therefore 
the main contribution to po comes from slowly varying 
terms in flee (i). 



To find the smooth terms flee (i) we write in Eq. ([5]) 
z„ = Yla/3 z a/3 (J n 13 ( an d similarly for z%) and notice that 
UH^^/Uit) oc exp(iu a pt). Here z afi =„(a\z n \P) n ; the 
matrix elements of are defined similarly. The inter- 
subband transition frequencies u) a p largely exceed the in- 
plane vibration frequencies, which are ~ to p and charac- 
terize time variation of the factors r~ 3 , in Hel\t). There- 
fore the first type of slowly varying matrix elements in 
Hel\t) are those containing only diagonal components 
z aa a" a and (z 2 ) aa cr" a . The second type comes from 
the terms oc z n z n i with n ^ n'. The corresponding 
slowly varying terms are proportional to \z a p\ 2 a^a^? 
with a ^ p. 

We will analyze the interaction H^) in the mean-field 
approximation. This approximation is justified when the 
number of nearest neighbors for each electron is large, 
so that correlations between out-of-plane states of neigh- 
boring electrons and their positions are averaged out. In 
other words, tracing over the positions of coupled elec- 
trons and over their out-of-plane states is done inde- 
pendently. In a strongly correlated 2D electron liquid 
the number of nearest neighbors is 6, on average, which 
makes the approximation reasonable. For typical times 
on the order of the duration of a collision with helium ex- 
citations or a radiation-induced inter-subband transition 
the relative shift of the positions of neighboring electrons 
is small. It is a good approximation then to describe the 
distribution of inter-electron distances r nn ' by the static 
pair correlation function g(r nn i). 

The many-electron frequency shift Aw^ of the 
\a) — > \(3) transition is determined by the ratio 

I Pfia- From the above arguments 



Tr 



?Y r~ 3 ,a a 7 p W 



} ~Fn 3 J 2 pp aPp , a ,, (10) 



where the prime over the sum indicates that n' ^ n and 



F = n. 



-3/2 



(11) 



Numerically, in a strongly correlated liquid F w 8.91—. 
From Eqs. and (fit))) with account taken 

of the commutation relation cr"^ , a", = 

$nn> (v%P'Sp',a> - (7^5 a ^\ the frequency shift of 

the |1) — > 1 2) transition due to the electron-electron 
interaction is 



Aw 2 i : 



Fe/n 



2 3/ 2 r 



2h 

ZoLQ.Pa.OL 



(^22 



2|zi2| 2 (pn - P22) 



2(z u 



Z22) 



(12) 



From Eq. (|12p. the frequency shift AW21 is determined 
by the populations of all subbands p aa - In Eq. (|12[) we 
dropped terms proportional to p\ a p a 2 with a ^ 1,2. 



These terms would be important if radiation excited 
transitions other than |1) — >• |2). However, such transi- 
tions are nonresonant, and for a comparatively small mi- 
crowave power that we consider they can be disregarded. 

The effect of the electron correlations on the position 
of the absorption line was considered earlier by Lambert 
and Richards^. Their results referred to the weak-power 
limit, p aa oc 5 a ,i, and the expression they were using 
differed from Eq. (|12[k in particular, the last term was 
missing all together. 

The physics of the terms proportional to the diago- 
nal and off-diagonal matrix elements of z, z 2 in Eq. (|12[) 
is different. The diagonal terms come from the fact 
that SEs in different out-of-planc states \a) have differ- 
ent static dipole moments in the out-of-plane direction. 
These static moments lead to the static Stark shift of 
the transition frequency u>2i of other electrons, which de- 
pends on the occupancies p aa . The term proportional to 
Z12P in Eq. (|12[) is dynamic in nature: it describes res- 
onant excitation transfer between different electrons^ 



D. Kinetic equation 

An equation for the single-electron density matrix for 
out-of-plane motion p a p can be obtained from the full ki- 
netic equation for po in the interaction representation^ 

Paa = (ipaPPP ~ laPPaa) + (S a ,l - ^afi) Impi 2 , 

pYl = -(iSu + 7o)/9l2 + -£fifl(/322 - /Oil), (13) 

where 

Soj = ujp — ix>2\ — Awgi- (14) 

In Eq. (fTB")) we do not incorporate explicitly the electron- 
electron interaction, which leads to establishing thermal 
distribution over lateral motion, with effective tempera- 
ture T e . We also do not incorporate the relaxation pro- 
cesses that lead to energy exchange between the electrons 
and the helium excitations. 

Parameters "/ a p in Eq. (|13|) are the rates of inter- 
subband transitions \a) \f3), which are essentially 
quasi-elastic. The electron momentum goes to helium 
excitations whereas the energy of motion normal to the 
helium surface goes into the kinetic energy of lateral mo- 
tion and the energy of the electron-electron interaction, 

/oc 
dtexp[i(e a - ep)t/H] 

x (exp[iqr n (t)] exp[-tqr n (0)]) . (15) 

Here, the overline indicates averaging over the thermal 
distribution (with temperature T) of the relevant excita- 
tions in helium; (...) denotes averaging over the in-plane 
electron motion with temperature T e . Expression (|15[) 
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FIG. 1: (color online) Sketch of the energy subbands and the 
electron transitions. Vertical (red) arrows indicate photon- 
induced transitions, while horizontal (blue) arrows indicate 
transitions due to elastic scattering by excitations in helium. 



is independent of the electron number n. In deriving it 
we used the aforementioned picture in which electrons 
are scattered by a short-range potential independently, 
that is one electron collides with a short-range scatterer 
at a time. However, during a collision the electron is in 
the in-plane potential created by other electrons, which 
is determined by ■ This potential is smooth on the 
electron wavelength and for ksT e 3> hu> p weakly af- 
fects short-range scattering in the absence of a magnetic 
fields 

The electron-electron interaction plays a critical role 
if the electron system is in a strong magnetic field B 
normal to the surface. Here, in the single-electron ap- 
proximation the electron energy spectrum is discrete (the 
Landau quantization), and the Born approximation used 
to derive Eq. (fTS")) would not apply. However, because 
an electron is in an a fluctuational electric field of other 
electrons Ef, its potential (and thus kinetic) energy is 
uncertain by eEfXx ~ hijj v , and if this uncertainty ex- 
ceeds huj c , where w c = eB/mc is the cyclotron frequency, 
the effect of the Landau quantization on scattering is es- 
sentially eliminated^ This condition was met in the ex- 
periment, even though the magnetic field was classically 
strong. 

The processes described by Eq. (fT5"j) arc shown 
schematically in Fig. [TJ The calculation of the rates 7 Q( g 
is similar to that of the in-plane relaxation rate in a cor- 
related electron system^ Even where T e ^ T, since the 
scattering described by Eq. (fT5j) is elastic and the elec- 
tron distribution is thermal, the transition rates satisfy 
the detailed balance condition 

lap = 7/3aCxp[(e Q - ep)/k B T e }. (16) 

The decay rate 70 of the off-diagonal matrix element 
P12 in Eq. (|13[) in the single-electron approximation for 
T e = T <C (e2 — ei)H was considered in Ref. |U Using the 
approach of Ref. [39| the results can be extended to the 
many-electron case and made applicable in the presence 
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of a magnetic field normal to the helium surface, 

70= i(5>ia + $>2a) + ^IVl ~ Vq 22 | 2 

/oo 
d£ (exp[iqr„(t)] exp[-iqr„(0)]) .(17) 
-oo 

The frequency detuning (|14[) depends on the distribu- 
tion of SEs, cf. Eq. (|12j) . which makes the overall quan- 
tum kinetic equation (|13[) nonlinear. As shown below, 
this nonlincarity may result in the bistability of the re- 
sponse of the electron system to resonant radiation. 

E. Absorption bleaching 

To find the stationary distribution of the system, 
Eq. (p~3|) must be complemented with the equation that 
describes electron energy relaxation. However, an im- 
portant conclusion about the distribution can be reached 
even without the analysis of energy relaxation in the case 
of weak to moderately strong microwave radiation, where 

« 70721- (18) 

Condition (fTSJ) is sufficient to suppress the conven- 
tional absorption saturation. Such saturation is well- 
known for two levels systems and requires that radiation 
makes the populations of the excited and ground states 
P22 and pn close to each other, so that the probabilities 
of radiation-induced transitions up and down in energy 
become close, too. In the saturation regime the popu- 
lation ratio pill P\\ significantly differs from the Boltz- 
mann factor exp[(ei — E2)/k^,T e ]. In our case instead of 
the levels we have subbands of electron motion. It is 
seen from Eq. (fT3|) that in the range (p~8|) the deviation 
of P21I p\\ from the Boltzmann factor is small even at 
resonance, Su> = 0, and even if one disregards transitions 
between subbands 1, 2 and other subbands. This shows 
that absorption saturation does not occur in our system. 

From Eq. (|13l) and from the detailed balance condition 
(Ti7J|) it follows that in the range (TT5|) the overall stationary 
distribution of the electron system both over the energy 
of lateral motion and over the subbands of out-of-plane 
motion is characterized by the same temperature T e , 

p aa MZ- 1 e-*-' h » T ; Z = Y J ?- ta/kBT °- (19) 

a 

Even though there is no absorption saturation in the 
range (TT5|) , the absorption significantly changes with the 
increasing radiation intensity. The mechanism of this 
change is absorption bleaching j 12 i 41 It can be understood 
from Fig. [TJ As a result of inter-subband scattering an 
electron that resonantly absorbed a photon in subband 
1 and made a transition to subband 2 can go back to 
subband 1, but with higher energy. It can then again 
absorb a photon. Now it can be elastically scattered from 
subband 2 to subband 3. It can also go back to subband 



1, again resonantly absorb a photon, and ne scattered 
into subbands 4, 5, etc. 

The cascade of photon-induced transitions and inter- 
subband scattering leads to population of higher and 
higher subbands and to a decrease of the population of 
subband 1. Such decrease, along with the fact that the 
populations of subbands 1 and 2 become closer with in- 
creasing T e , lead to an overall decrease of resonant ab- 
sorption. 

F. Bistability of resonant response 

A complete analysis of nonlinear response requires in- 
corporating inelastic processes that lead to energy ex- 
change between the electrons and excitations in helium. 
By now several types of inelastic processes have been 
identified, including the weak inelasticity of scattering by 
helium vapor atoms, the inelastic two-ripplon scattering 
and the inelastic scattering by phonons in liquid helium. 
For all of them, the energy relaxation rate depends on 
the electron state and helium temperature^ i 14 ' 16 i 43 

The description of energy relaxation is simplified by 
the fact that it is slow. Between inelastic scattering 
events electrons have time to exchange energy with each 
other, and the state of the electron system is described 
by electron temperature. Therefore one can introduce a 
single energy relaxation rate i/g, which is independent 
of the electron state and depends only on T e ,T, and 
— — — . This leads to an implicit equation for T e . 
It has the form of a simple energy balance equation 

v E k B {T e — T) = hjj F rZ- 1 [e- tl,k ^ - e" £2/fcBTc ], 
r = 7j^7o/[7o 2 + ~ "ai - A ^2i) 2 ]- (20) 

We have introduced here absorption coefficient r, which 
is determined by Eqs. (fT3| and (|14|) . The left- and right- 
hand sides of Eq. ([20|) are, respectively, the energies dis- 
sipated into helium and absorbed by an electron from 
radiation per unit time. For all energy relaxation mech- 
anisms considered so far, the left-hand side of Eq. (|20|) 
is an increasing function of T e — T, at least for not too 
large T e -T. 

Parameter r is a nonlinear function of T e , as seen from 
Eqs. {H]) and (fT3|). Therefore Eq. (J5DJ) is fairly compli- 
cated. A convenient way of understanding the possible 
types of its solution is based on the graphical solution, 
an example of which is shown in Fig. [5J The plot shows 
the rates of energy absorption and relaxation calculated 
for T = 0.4 K where the scattering is dominated by in- 
teraction with helium vapor atoms. In this regime, the 
T e -dcpcndcncc of ve, 70 and AW21 are calculated for the 
typical experimental parameters with account taken of 
200 subbands of electron motion normal to the helium 
surface. 

The left-hand side of Eq. (f2TJ)) is shown in Fig. [5] by a 
dashed line, while the right-hand side (rhs) is shown by 
the solid lines for ujp— = 0.1 and 0.4 GHz (lines a and 
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b, respectively). The major part of T e -depcndence in the 
rhs comes from the frequency shift Aw 2 i- From Eq. ([T^|) . 
Aw2i rather quickly increases with the increasing T e . In 
a certain range of up — U21, the change of T e can bring 
AW21 in resonance with u>p—u>2i- This leads to a resonant 
peak of r and thus of the rhs of Eq. (|20|) as function of 
T e seen in Fig . gj 

The heights of the resonant peaks in Fig. [5] depend on 
up — oj-21- From Eq. the coefficient in front of r 

decreases with increasing T e . The maximal value of r at 
resonance, r max oc I/70, also decreases. Therefore the 
higher is the value of T e where AW21 = — U21, the 
smaller is the height of the peak of the rhs of Eq. ([20)) . 
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FIG. 2: (color online) Graphical solution of the energy bal- 
ance equation ([20jl for T = 0.4 K and 7i a = 4.0 x 10 7 cm 2 . 
The dashed line is the energy loss rate [the left hand side of 
Eq. (|20p ] due to inelastic scattering by helium vapor atoms. 
The solid lines show the rate of the microwave power ab- 
sorption [the right hand side of Eq. (|20p ] calculated for 
ujf = 104.5 GHz, Q.r = 2.0 MHz, and u F - oj 2 i = 0.1 GHz 
(a) and 0.4 GHz (b). For the case (a), the stable solutions 
of Eq. (|20p . which are given by the intersections of solid and 
dashed lines, are marked by (red) circles. 

The solutions of Eq. (f2T)| are determined by the cross- 
ings of the dashed and solid line in Fig. [2] They give the 
effective temperature T e in the stationary state. Where 
the microwave frequency is relatively far from resonance, 
i.e., top — W21 is relatively large, there exists only one so- 
lution for T e , with T e close to the helium temperature T. 
Closer to resonance, on the other hand, Eq. ([2U|) can have 
two stable solutions with comparatively large and small 
T e —T. They are marked by circles in Fig. [2] In this case 
Eq. (f2"0"|) also has a stationary solution for intermediate 
T e , as seen from Fig. [5] By writing the energy change 
per unit time as the difference of the right- and left-hand 
sides of Eq. ([20)) . one can see that this intermediate so- 
lution is unstable. 

The left and right hand sides of Eq. (|2T))) display a 
qualitatively similar behavior for other parameters and 
for other relaxation mechanisms. This indicates that the 



onset of bistability is a general consequence of absorption 
bleaching. 

The occurrence of two stable regimes can be under- 
stood in the following way. In one regime, which cor- 
responds to lower T e , the overall detuning of the elec- 
tron transition frequency U21 + AW21 from the radiation 
frequency uf is large compared to 70. Therefore, the 
absorption rate is comparatively small and T e is self- 
consistcntly low. On the other hand, if T e is sufficiently 
high, W21 + AW21 is close to wj?, absorption is compar- 
atively strong and T e is self-consistently high. We note 
that, since AU21 is positive, such bistability can occur 
only for cup — U21 > 0. 

The bistability leads to the hysteresis of electron ab- 
sorption with varying up — U2\. If we start from large de- 
tuning ujp —UJ21 (curve A in Fig. [2]) and then decrease it, 
the system will stay in the small-absorption state, which 
corresponds to left circle in Fig. [2j On the other hand, 
if we start from small up — U21, the system will be on 
the decreasing with increasing T e part of the resonant 
solid line in Fig. [2] As up — U21 increases, the system 
will get into the region of strong absorption, but ulti- 
mately will switch to the low-absorption branch for large 
up — U2i- Below we present the experimental observation 
of frequency shift, bistability, and hysteresis and make a 
comparison with the theory. 



III. EXPERIMENT 

The details of the experimental apparatus have been 
described elsewhere^. SEs were accumulated on the sur- 
face of liquid 3 He placed about half way between two 
round electrodes of diameter 26 mm separated by ap- 
proximately 2.6 mm and forming a parallel-plate capac- 
itor. A positive voltage Vb was applied to the bottom 
electrode to create an electric field E± perpendicular to 
the helium surface. In addition, a magnetic field B was 
applied perpendicular to the surface using a supercon- 
ducting magnet placed around the experimental cell. Mi- 
crowave radiation at fixed frequency uf/2it = 104.5 GHz 
was passed through the experimental cell, and the inter- 
subband splitting 62 — 61 was brought in resonance with 
hup by varying E±. 

The top plate of the capacitor consisted of two 
concentric-ring electrodes (the Corbino disk geometry), 
which were used to measure the longitudinal magneto- 
conductivity of the electron system a xx . The Corbino sig- 
nal proportional to o xx was recorded while slowly vary- 
ing Vb at a fixed microwave power and at a fixed value of 
the magnetic field B. The variation of u xx caused by the 
radiation-induced heating of SEs allowed us to observe 
resonant response from inter-subband absorption. Sub- 
tracting the value of the signal recorded in the absence 
of radiation, the relative change Aa xx due to the 

heating was determined and plotted vs Vb. The mea- 
surements reported here were done in the temperature 
range from 0.2 to 0.4 K. 
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For each curve Aa~^/a~^ vs Vb, the electron tem- 
perature T e was calculated from the previously discussed 
dependence on T e of the quasi-elastic electron transport 
time Ti 10 ' 14 For magnetic fields B used in the experiment, 
the many-electron conductivity a xx follows the Drude 
law^ 



previously ; 10 ' 14 in this case radiation-induced heating 
leads to a decrease of the relaxation time r. Therefore 
a~} should decrease at resonance. 



(l + fx 2 B 2 )/a , 



(21) 



where Co = n s efj. is the conductivity in zero magnetic 
field and \i = er/m is the electron mobility. For [iB 1, 
a xx is proportional to r. To find the proportional- 
ity coefficient, which depends on the electron density 
n s , the Corbino signal was recorded without radiation 
while SEs were slowly cooled down until they formed a 
Wigner crystal. The crystallization could be easily de- 
tected as an abrupt change in a xx £^ The value of n s 
was found from the crystallization temperature and the 
well established critical value of the plasma parameter 
r = e 2 {-Kn s f/ 2 /k B T^ 130. 

The experimental procedure for direct absorption mea- 
surements was similar to that described in Ref. 0. The 
power of the microwave radiation that passed through the 
cell was measured with an InSb bolometer mounted in- 
side the cryostat. Resonant absorption could be observed 
as the variation of the bolometer signal at T = 0.4 K as 
electrons were driven through resonance by sweeping the 
bottom plate voltage Vb- Simultaneously, this voltage 
was modulated with a small sinusoidal signal at frequency 
10 kHz, and the in-phase and quadrature components of 
the demodulated bolometer signal were obtained using a 
dual-phase lock-in amplifier. With such an arrangement, 
in the limit of small modulation amplitude the in-phase 
component of the obtained signal is proportional to the 
derivative of the power absorbed by SEs with respect to 
Vb- Normally the quadrature component of the signal 
should be zero. The absorption lincshapc was obtained 
by numerical integration of the signal with respect to Vb . 



IV. RESULTS 
A. Resistivity measurements 

In Fig. [3] we show the relative radiation-induced con- 
ductivity change — Aa'^ / & xx as a function of the voltage 
that presses the electrons to the helium surface Vb- All 
plots were obtained from the data taken at T = 0.4 K and 
electron density n s , = 4 x 10 7 cm~ 2 . The radiation power 
was measured at the output of the microwave source; 
shown in Fig. [3] is the ratio of the applied power (in deci- 
bel) to the maximum power, which was approximately 
1 mW. Experimental traces were obtained by increasing 
Vb slowly to drive SEs through resonance. A conversion 
of the increase of Vb to the increase of the transition 
frequency w 2 i by « 25 MHz is discussed in Sec. V A. 

At T = 0.4 K, electron scattering is mostly determined 
by collisions with helium vapor atoms. As we showed 
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FIG. 3: Relative change of the reciprocal conductivity 
—Aa~x/<Jxx vs the voltage applied to the bottom electrode 
Vb obtained for SEs on liquid 3 He at T = 0.4 K, n 3 = 
4 x 10 7 cm" 2 and several values of the input power of ra- 
diation. The lines a to k correspond to -17.5, -14.0, -11.4, 
-8.8, -6.2, -4.2, -2.3 and -0.6 dB of the input power expressed 
as a fraction (in decibel) of the maximum power ~ 1 mW. 



The resistivity curves in Fig. [3] show a pronounced res- 
onant peak. The height of the peak monotonically in- 
creases with the input power. At small power levels, the 
resistivity shows a slight increase on the sides of the res- 
onance (respectively, — Aa~^ becomes negative). This 
effect is attributed to the contribution of ripplon scat- 
tering to the momentum relaxation r&teJ^ As the power 
increases, the resonance shifts toward lower values of Vb, 
which corresponds to the increase in the transition fre- 
quency. Also, the lineshape becomes asymmetric. We 
believe the shift is due to the many-electron effect de- 
scribed by Eq. (fl"2"|) and is associated with thermal pop- 
ulation of the excited subbands as the effective electron 
temperature increases. 

At high power levels, when the voltage is swept 
through the resonance from the 1ow-Vb side, the re- 
sistivity jumps abruptly to a lower value; respectively, 
—Aa xx /a~ x jumps up. This unusual behavior is demon- 
strated in Fig. |4j The data used for this plot were taken 
at similar conditions to those in Fig. [3J except that the 
sweeping rate was slower by a factor of two. The jump 
always occurs on the low-field side of the resonance and 
is well reproducible. When the sweeping direction is 
reversed, the resistivity curve first follows the forward- 
sweeping curve as it passes through the resonance. How- 
ever, on the low-field side it does not show any jump but 
slowly goes to zero with decreasing Vb- This behavior 
is also well reproducible. Therefore, overall the system 
displays a well-pronounced hysteresis. 

The resistivity jump becomes more prominent and the 
width of the hysteresis loop increases at lower tempera- 
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FIG. 4: Relative change of the reciprocal conductivity 
—Aa^x /o~xx vs Vb obtained at T = 0.4 K, n s = 4 x 10 7 cm" 2 
and at the maximum input power (OdB). The dashed and solid 
lines correspond to sweeping through the resonance from the 
low and high side of Vb, respectively, as indicated by the ar- 
rows. 

tures. In Fig. [5] we plot Aa~^/a~^ vs Vb for the data 
taken at T = 0.2 K, n s = 4.2 x 10 7 cm" 2 , and at three 
different input power levels. At this temperature, the 
relaxation rate is determined by the interaction of the 
electrons with ripplons, and electron heating leads to an 
increase of the relaxation time r. Therefore a~% increases 
at resonance^ In addition to the abrupt jump and hys- 
teresis observed at T = 0.4 K, the lineshape shows an 
interesting structure. It seems to be a superposition of 
a broad peak centered at Vb ~ 11-9 V and a high peak 
located on the 1ow-Vb side of the broad peak. While this 
latter peak grows and shifts with the increasing input 
power, the position of the broader peak does not change 
and the height of the peak slightly decreases with the in- 
creasing power. We identify the high peak as being due 
to resonant transitions from the ground state to the hrst 
excited state. The broader peak remains unexplained at 
this time. 

In addition to the peaks described above, a small satel- 
lite peak appears on the low-field side of the resonance. 
This peak becomes visible only at very high power lev- 
els. It was reported in our previous worb^ and was at- 
tributed to the excitation of collective in-plane electron 
vibrations 4^ At extremely high radiation intensities, we 
observe a number of satellite peaks on the both sides 
of the resonance. Some of them can be due to reso- 
nant transitions between higher excited subbands and 
emerge through self-sustained absorption, which was re- 
cently discovered in our experiments.— The origin of 
other peaks is not clear at the moment. A detailed de- 
scription of the satellite structure of the resonance will 
be given elsewhere. 

The signal strongly changes at lower electron den- 
sity. In Fig. [6] we plot Aa~lja~l vs Vb for the data 
taken at the same temperature as in Fig. [5] but at n s = 
8.0 x 10 6 cm" 2 . Here, the intensity (area) of the broad 
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FIG. 5: Relative change of the reciprocal conductivity 
Aa^/a'^ vs V B obtained at T = 0.2 K, n s = 4.2 x 10 7 cm" 2 
and three different values of the input power: -6.2, -5.4 and 
-4.2 dB (lines a, b and c, respectively). The dashed and solid 
lines correspond to the sweeping through the resonance from 
the low- and high-VB sides, respectively. A satellite peak that 
appears on the low-field side of the resonance is indicated by 
the full triangle. 
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FIG. 6: Acr^/o-^ 1 vs V B for T = 0.2 K, n s = 8.0 x 10 6 cm" 2 , 
and -2.3 dB of the input power. The dashed and solid lines 
correspond to the sweeping through the resonance from the 
low and high Vb side, respectively, as indicated by arrows. 
Inset: plots obtained under the same conditions, but with 
the slower by a factor of two rate of sweeping from the lower 
(opened squares) and the higher (solid squares) Vb side. 



peak is significantly smaller, while the higher peak is 
much narrower. Interestingly, for lower electron density 
the signal shows abrupt jumps for both sweeping direc- 
tions. It has a characteristic sawtooth shape expected for 
switching between two stable states. The abrupt switch- 
ing for sweeping the voltage up and down is highlighted 
in the inset of Fig. [51 



10 



B. Absorption measurements 

We have also studied absorption directly by measur- 
ing the power of the microwave radiation passed through 
the electron system. The experiment was done by si- 
nusoidally modulating the potential Vb and sweeping 
the central point Vb about which the modulation was 
performed. The measured transmitted power was then 
Fourier-transformed using a lock-in amplifier. Unexpect- 
edly we found that, in the absorbed power, along with 
the standard in-phase component /(Vb), for high radia- 
tion power the signal displayed a quadrature component 
Q(Vb)- The data were obtained for the modulation fre- 
quency w m /27r=10 kHz and the modulation amplitude 
V m = 10 mV, which corresponds to the modulation am- 
plitude pa 25 MHz of the transition frequency ui 2 i (see 
below). For modulation oc cos u) m t we define, following 
the standard notations, the quadrature component as a 
coefficient at — sinoj m t. 

In Fig. [7J we show functions 



/Vb r 
dViHVi), Q(Vb) = J dV B Q(V£), 

(22) 

sometimes called integrated lineshapes. They were ob- 
tained by sweeping V B upward for T = 0.4 K and 
n s = 4 x 10 7 cm" 2 at different power levels. The lower 
limits of the integrals in Eq. ([22]) were chosen at a value of 
Vb well below the resonance; the result weakly depended 
on this limit, as it is clear from Fig. [JJ For T = 0.4 K the 
intrinsic lincwidth of the absorption spectrum is rather 
small and the observed lineshapc is determined by the in- 
homogeneous broadening due to the nonuniformity of the 
electric field at the helium surface as well as the nonuni- 
formity of the microwave field 7 -^. 

As expected, for low radiation power the quadrature 
component of the absorption signal is zero. However, at 
higher powers hysteresis occurs within each modulation 
cycle itself. This gives a distinctive quadrature compo- 
nent in the modulated output, sec Sec. V. The effective 
lineshape which is obtained from the integrated modu- 
lation signal then has a quadrature component and an 
offset in both components. Note that the quadrature 
component itself is not hysteretic - it is the same sweep- 
ing the voltage up or down. But this component gives a 
sensitive indication of local hysteresis in each modulation 
cycle. The features of the integrated lineshape become 
more prominent with the increasing power. For lower T, 
the quadrature component is seen for lower power levels. 



V. DISCUSSION 

A. The many-electron frequency shift 

The resonance condition for exciting an electron from 
the ground to the first excited level of motion normal to 
the helium surface is determined by the many-electron 
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FIG. 7: Integrated lineshape of the absorption signal mea- 
sured with SEs on liquid 3 He for T = 0.4 K, n s = 4x 10 7 cm" 2 
at different power levels. The graph shows the in-phase (panel 
A) and the quadrature (panel B) components of the measured 
signal. The lines a to e correspond to -27.7, -20.6, -17.5, -14.4, 
and -11.4 dB of the input power attenuation. 



frequency shift, uip — U21 = AW21 with Aw2i(T e ) given 
by Eq. [12] The position of the resonance depends on 
the electron temperature T e . In turn, T e depends on the 
absorbed power and for fixed incident power is expected 
to be maximal at resonance. In the parameter ranges 
used in our experiment <j~£ varies monotonically with T e . 
Therefore we assume that the maxima of the resistivity 
curves plotted in Fig. [3] correspond to exact resonance. 
Then, for each power level, the frequency shift AW21 can 
be found from the shift of the maximum of the corre- 
sponding curve with respect to the resonant absorption 
curve measured at very low power using the InSb bolome- 
ter. For this curve one can assume T e = T, and then in 
Eq. (|12p for AW21 one can set p aa = <5 Qj i. Generally, the 
parameters z a p in Eq. (|12|) depend on the pressing field 
E± and thus on Vb; however, this dependence is smooth, 
and in the narrow range of Vb that we studied it could 
be disregarded. 

To find the shift in frequency units, the frequency of 
the microwave source was changed by a known amount 
and the corresponding shift in Vb was recorded. Such 
measurements were done at low power so that the many- 
electron shift remained constant. This procedure allowed 
us to establish the relationship between U21 and Vb- The 
conversion factor, i.e., the slope of W2i/27r vs Vb, was 
found to be 2.5 ±0.1 GHz/V. 

In Fig.[8]we plot the many-electron frequency shift and 
the electron temperature determined for the maximum of 
each curve in Fig. [3] We remind that T e was obtained by 
comparing the measured with its theoretical value 
for given n s and T. The frequency shift is counted off 
from its value for T e = T . It monotonically increases 
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FIG. 8: Many-electron frequency shift Acj2i vs electron tem- 
perature T e obtained from the data shown in Fig. [3] (squares); 
Aoj2i is counted off from its value for T e — T = 0.4 K . The 
solid line shows the result of calculations using Eq. [12] for 
n s — 4 x 10 7 cm -2 . The dashed line shows the calculation 
that excludes the term proportional to |zi 2 | 2 in Eq (fl2|) . 

with T e and reaches the value of about 1 GHz for the 
highest-power curve in Fig. [3] For comparison, we also 
show in Fig. |S] the frequency shift calculated from Eq.[T2l 
for n e = 4 x 10 7 cm~ 2 . In this calculation we limited 
the number of the states of out of plane motion to 20, 
which was sufficient for the studied range of T e . At high 
T e , the shift observed in the experiment is about 20% 
larger than our theoretical estimate. This discrepancy 
might be attributed to the nonuniformity of the electron 
density and the microwave power across the electron layer 
and also to the corrections to the approximations used to 
obtain Eq. [121 and to evaluate T e . 



B. Nonlinear absorption and bistability 

The experimental results on the bistability and hys- 
teresis of nonlinear response discussed in Sec. IV (sec 
Figs. [Uini) are in qualitative agreement with the theory 
presented in Sec. II. The data clearly demonstrate hys- 
teresis: the state of the system depends on the direction 
of sweeping Vb or, cquivalently, the single-electron tran- 
sition frequency u>2i- The abrupt jump in Aa~^/cr~^ 
observed in the experiment is due to switching between 
different states of the electron system that coexist in a 
limited parameter range. Switching occurs once the cor- 
responding state disappears with the varying parameter. 

Sharp jumps in Figs. Hl-JfJ] were observed with increas- 
ing W21 for switching from the states with comparatively 
small absorption and low T e to high-T e states. Such 
switching is expected from the theory and can be inferred 
from Fig. [2] If one thinks of moving curve B in Fig. [5] to 
the left, one can see that the low-T e state will disappear. 
For data taken at T = 0.4 K (see Fig. [J}, T e jumps from 



about 0.4 K to about 3.1 K. For data taken at T = 0.2 K 
and with similar electron density (see Fig. [5]), the jump 
is from about 0.3 K to about 2.6 K. 

For the reversed sweeping direction, from higher to 
lower UJ21, the response displays an abrupt jump only 
for low electron densities, see Fig. [TJ] For higher densities 
a sharp switching to the small- absorption state was not 
observed, instead the response smoothly varied with the 
parameters in the range of hysteresis. This indicates that 
the dynamics of the system is more complicated than the 
model of Sec. II suggests. In particular, the spatial in- 
homogeneity of the electron density and the microwave 
field should play a role. 

The inhomogencity would favor a sharper upward 
jump of the absorption. One can think that, with the 
increasing Vb, there is formed a "critical nucleus" of 
high absorption and high T e , which then quickly expands 
throughout the system. On the other hand, with the de- 
creasing Vb , the areas that exist only in the low T e -state 
gradually increase in size. Regrettably, these processes 
are hard to characterize quantitatively without a detailed 
picture of the radiation power distribution. Another in- 
dication of a more complicated nature of the system is 
the broad peak of Aa~^ /a~^ vs Vb in Figs. [5] and 
which the model of Sec. II does not explain. We note 
that the narrowing of the hysteresis loop for lower den- 
sity in Fig. [TJ] compared to Fig. [5] can be understood by 
noticing that AW21 scales with the electron density as 
n 3 / 2 for the same T e . This gives a factor of 12 difference 
between the values of Au>2i for the densities in Figs. [5] 
and[?5] 



C. Absorption signal for modulation across the 
hysteresis loop 

The theory of Sec. II also explains the complicated 
behavior of the microwave power absorption described in 
Sec. IV. B, see Fig. [7J This type of behavior was also 
reported by Glasson et al. for SEs on liquid 4 He and was 
attributed to hysteresis of microwave absorption* 2 ^ In the 
absorption experiment, the transition frequency W2i(i) is 
swept up and down with a modulation period To. Within 
a modulation period, the absorption for increasing and 
decreasing o->2i(i) can be different. This happens if the 
modulation amplitude is large enough, so that ui2i{t) goes 
across the whole region where the two stable states of the 
electron system coexist, i.e., across the hysteresis region. 

We now show that hysteresis leads to the behavior 
seen in Fig. [7] Suppose the frequency is modulated si- 
nusoidally, W2i(t) — <I>2i = Acos(uj m t), where the overlinc 
means period average. In the experiment, W21 is deter- 
mined by the value of the potential Vb about which the 
modulation is performed, W21 = ^2i(Vb). For the chosen 
modulation phase, uj2i{t) is swept down during the first 
half of the modulation period 7r/w m and up during the 
second half. 

If W2i(£) goes across the hysteresis region, the absorbed 
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microwave power S(tJ2i(t)) takes on different vaiues, 
S-(uj2i(t)) and S+(w2i(i)) on the down and up sweeps, 
rcspectiveiy; note that functions S± have discontinuities 
where the absorption switches between the high- and low- 
absorption branches, which we denote by <S>(w 2 i) and 
S<(u>2i), respectively. Where U)2i(t) goes across the hys- 
teresis region, function S— first evolves along the branch 
5> and then switches to S< , whereas function S+ evolves 
along S < and then switches to 5> . 

The in-phase and quadrature components of the ab- 
sorbed power detected with a lock- in amplifier, I = /(Vb) 
and Q = Q(Vb), are 

Tr/uj m 

1=^- J COs{0J m t)[S-{U2l{t))+S + {U2l{t))]dt, 


7r/u; m 

Q = ^ / M"mt) [S+(uj 21 (t)) - S-Mt))] dt.(23) 



From Eq. ([23]) . where the modulation cycle goes across 
the hysteresis loop the signal has a nonzero quadrature 
component, 

i r A 

Q = —r [S+(lD 2 i + x) - S-(Q 2 i+x)]dx. (24) 

The value of Q does not depend on what part of the 
modulation cycle is covered beyond the hysteresis region, 
as outside this region S-\~ — S— . 

We now explain the behavior of the experimentally 
measured parameters X and Q in Fig. [7] where, as Q\2 is 
increased, it goes through the hysteresis region in such a 
way that the system switches from the lower- to a higher- 
absorption branch. We start with the quadrature com- 
ponent Q. Let us assume that hysteresis occurs in the 
region lol < u>2i < ujj. If the modulation amplitude 
A < (ujh ~ujl)/2 then as the voltage Vb is swept up, the 
absorption will instantaneously switch from the low to 
the high branch. In this case the quadrature component 
Q = 0. But if A > (ojh — wl)/2, the hysteresis loop will 
occur in each modulation cycle, in the appropriate range 
of Q 2 i (%)■ Then 



Q = Q(Vb) = (irAy 



dx[S < (x) — S > (x)] 



x [0(a>2i -ujh + A)- 0(u>2i A)} , (25) 

where 0{x) is the Hcaviside step function and Q21 = 
W2i(Vb); we used that S+(ui) = S < (u>) and S-(ui) — 
S>(u>) for u>l < oj < ujh- 

From Eq. (f2"5|) . Q(Vb) remains constant for ujh — A < 
u>2i < w l + A. This constant is negative and is equal to 
—Sn/irA, where S\j is the area of the hysteresis loop on 
the plane (u2i,S). Respectively, the integrated quadra- 
ture Q(Vb) linearly decreases with increasing Vb in this 
range and becomes constant once 0)21 (Vb) reaches lul+A, 



in agreement with Fig. [7] This constant, which is the off- 
set of the integrated quadrature AQ on the high-Ve side, 
is given by 



AQ = -(2CV/7r)S n [1 - (ujh - lo l )/2A] , 



(26) 



where Cy = (du>2i/dVB) 1 is the reciprocal slope of tran- 
sition frequency vs voltage. 

The offset of the integrated quadrature is a clear in- 
dication of hysteresis in the system, which thus can be 
revealed by studying the response for a finite modulation 
amplitude. The increase of the slope of the curves in 
Fig. [7] with increasing power corresponds to the increase 
of the area of the hysteresis loop S\j ■ 

We now consider the in-phase component. If function 
S(ui) were single-valued and smooth, S(u>) = S-(u>) — 
S+(ui), in the limit of small modulation amplitude A — > 
we would have /(Vb) ~ Ad u S, which gives I(Vb) « 
ACy S (u)2i(Vb)) ■ Then X(Vb) goes to zero on the both 
sides of the absorption peak. 

The behavior of I becomes different in the pres- 
ence of the hysteresis loop. We assume that the di- 
rection of sweeping Vb is such that we move along the 
low-absorption branch and jump to the high-absorption 
branch, that is, in our system, we increase Vb- In 
the limit A — > (in particular, A <C <jJji — wj,), we 
have I(Vb) ~ ^4CVS<(w2i(Vv)) on the low-absorption 
branch 5<(w2i), i.e., for U21 < ojh- After the jump to 
the high-absorption branch 5>(o;2i) we have, as before, 
/(V B ) w A[d w S>] W2i . Then, for w 2 i > w H , X(V B ) w 
ACy [S > (ui2i(Vb)) - AS 1 ], where AS* = S>(w H ) - 
S<(u)h — 0) is the height of the jump of S at ujh- This 
shows that, as Vb will have gone over the resonant ab- 
sorption peak, where «S> (w 2 i(Vb)) — > 0, I will become 
negative. 

The offset AI on the high-Vs side is given by the 
height of the jump of the absorption coefficient, AX = 
— ACy AS . Such offset is a characteristic feature of the 
hysteresis. It sensitively depends on the interrelation be- 
tween the modulation amplitude and the width of the 
hysteresis loop^ 

We now discuss the case where the modulation range 
2 A exceeds the width of the hysteresis region and assume 
that this region is narrow. To the leading order in ui H — 
u>l and in A, for uil + A > u>2i > — A we have 
/(Vb) ~ 7t _1 AS[sinui m tL + sinw m i#], where the values 
of ti,,H are given by equations u)2i+Acosu m tL.H — wl,h 
with < u) m tL,H < 7T. This order of magnitude estimate 
is obtained by disregarding the change of <S> , S < within 
the hysteresis loop. 

The positive value of / reduces the negative offset of 
T. Integrating the above expression for /(Vb) over Vb in 
the range where uil + A > <D 2 i(Vb) > — A, we obtain 
for the offset 



AT w n^ACvAS 



K(l - H 2 ) 1/2 - arccosH 



H= 1 - A-\lj h -u l ). 



(27) 
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For a narrow hysteresis loop, <C 1, the offset 

is small, it scales as [(ujh — w l)M] 3 ' 2 - In the central part 
of the region < ujh—^l < 2 A the offset weakly depends 
on A and is almost linear in ujh — (jJl- As ujh ~ w l ap- 
proaches 2A the offset smoothly approaches the small- A 
value — ACyAS, which applies to the leading order in A 
for ujh — ijJl > 2 A. The reduction in offset compared to 
the small- A limit is easy to understand, since going round 
the hysteresis loop for 2A > ljjj — the system sam- 
ples both absorption branches in each modulation cycle. 
Both theoretically and experimentally, the overall neg- 
ative offset increases with increasing microwave power, 
as both the height of the absorption jump AS and the 
width of the hysteresis loop increase. 

VI. CONCLUSIONS 

This paper describes the theory and the experimen- 
tal observation of the long-sought intrinsic optical bista- 
bility and hysteresis in a quasi two-dimensional electron 
system. The effect is a result of strong electron corre- 
lations and the hierarchy of the relaxation times, where 
the electron-electron momentum and energy exchange is 
the fastest process, followed by a much slower momentum 
relaxation due to scattering by the random short-range 
potential of ripplons or helium vapor atoms, followed by 
still much slower energy relaxation due to electron en- 
ergy exchange with the environment. As a result of this 
hierarchy, a moderately strong resonant excitation of the 
transitions |1) — ^ |2) between the lowest subbands of 
motion normal to the surface leads to redistribution of 
the electrons over the subbands. This causes absorption 
bleaching due to depletion of the lowest electron subband. 
It occurs for much smaller power than that required for 
absorption saturation. 

Because of the electron-electron interaction, the redis- 
tribution of the electrons over the subbands leads to a 
change of the frequency UJ21 of the |1) — > 2) transition. 
This causes bistability of the irradiated system. For a 
given radiation frequency ujp the frequency LO21 can be 
"tuned" close to Wf as a result of relatively strong elec- 



tron population of excited subbands. Then absorption 
is indeed strong and the excited subband population is 
self-consistcntly significant. Alternatively, the resonance 
can be not that good, the electrons absorption is not that 
strong, the electrons mostly stay in the lowest subband, 
and W21 is self-consistently detuned from ujp. The de- 
veloped microscopic theory of the many-electron system 
provides a full account of this behavior. 

Experimentally, for electrons on 3 He we observed hys- 
teresis of microwave absorption. The observation was 
made by measuring the magnetoconductivity in the 
Corbino geometry. The data refer to the regime of strong 
electron correlations. Different temperatures were stud- 
ied, which made it possible to investigate the scattering 
primarily by the helium vapor atoms and by ripplons. 
The data for different radiation power show the increase 
of the hysteresis loop with the power. A quantitative 
agreement is obtained between the data and the theory 
of the electron-electron interaction induced change of W21 
as function of the electron temperature. 

We have also observed absorption hysteresis by di- 
rectly measuring the transmitted microwave power near 
the resonant frequency W21, which is modulated though 
the applied holding field. The hysteresis is manifested as 
a characteristic offset of the integrated lineshape in the 
in-phase component of the modulated response and by 
a distinctive quadrature component at finite modulation 
amplitude. 

The results of this paper provide a new insight into 
the role of the electron-electron interaction in strongly 
correlated systems and present new types of resonant 
nonlinear phenomena in quasi two-dimensional electron 
systems. 
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